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We report muon spin relaxation (pSR) measurements us- 
ing single crystals of oxygen-intercalated stage-4 La2CuO4.11 
(LCO:4.11) and Lai.ssSro.izCuCU (LSCO:0.12), in which neu- 
tron scattering studies have found incommensurate magnetic 
Bragg reflections. In both systems, zero-field fiSK measure- 
ments show muon spin precession below the Neel temperature 
Tjv with frequency 3.6 MHz at T — > 0, having a Bessel func- 
tion line shape, characteristic of spin-density-wave systems. 
The amplitude of the oscillating and relaxing signals of these 
systems is less than half the value expected for systems with 
static magnetic order in 100 % of the volume. Our results 
are consistent with a simulation of local fields for a heuristic 
model with: (a) incommensurate spin amplitude modulation 
with the maximum ordered Cu moment size of ~ 0.36 /j,b', (b) 
static Cu moments on the CUO2 planes forming "islands" hav- 
ing typical radius 15 ~ 30 A, comparable to the in-plane su- 
perconducting coherence length; and (c) the measured volume 
fraction of magnetic muon sites F M increasing progressively 
with decreasing temperature below Tjv towards F M ~ 40 % for 
LCO:4.11 and 18 % for LSCO:0.12 at T -» 0. These results 
may be compared with correlation lengths in excess of 600 A 
and a long range ordered moment of 0.15 ±0.05 fis measured 
with neutron scattering techniques. In this paper we discuss a 
model that reconciles these apparently contradictory results. 
In transverse magnetic field /iSR measurements, sensitive to 
the in-plane magnetic field penetration depth X a b, the results 
for LCO:4.11 and LSCO:0.12 follow correlations found for un- 
derdoped, overdoped and Zn-doped high-T c cuprate systems 
in a plot of T c versus the superconducting relaxation rate 
a(T — * 0). This indicates that the volume-integrated value 
of n s /m* (superconducting carrier density / effective mass) 
is a determining factor for T c , not only in high-T c cuprate 
systems without static magnetism, but also in the present 
systems where superconductivity co-exists with static spin- 



density-wave spin order. Q 



I. INTRODUCTION 

The interplay between superconductivity and mag- 
netism is one of the central issues of high-T c supercon- 
ductivity (HTSC), which has been extensively studied 
both experimentally and theoretically. [1] . In particular, 
dynamic [2-8] and static [9-16] spin correlations with in- 
commensurate wave vectors have been found by neutron 
scattering. These, as well as X-ray photoelectron studies 
[17-19], have been discussed in terms of a "stripe" mod- 
ulation of the spin and charge densities [20-23]. Yet it 
has not been clear whether static magnetism supports 
or competes with superconductivity, or if the magnetism 
and superconductivity coexist in the same microscopic 
regions of the Cu02 planes or in different regions. While 
neutron measurements make it clear that long range mag- 
netic order exists, the Bragg peak intensity integrates 
over the sample volume and is not, therefore, sensitive to 
microscopic spatial variations in the order parameter. 

Muon spin relaxation (/iSR) measurements [24,25] pro- 
vide a complementary probe in this regard. In a mag- 
netic material having a heterogenous structure, such as 
"magnetic" and "non-magnetic" regions, /iSR data are 
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composed of two different signals, corresponding to dif- 
ferent environments with signal amplitudes roughly pro- 
portional to their volume fractions. Local magnetic fields 
at muon sites result primarily from the dipolar interac- 
tion. In antifcrromagnctically ordered systems the local 
field decays very quickly with increasing distance from or- 
dered spins, resulting in an effective range of the field of 
<~ 10-f5 A in cuprate systems. Thus, any heterogeneous 
magnetic structure in HTSC, having a length scale larger 
than this, should produce multiple and distinguishable 
ytiSR signals. 

In this paper, we report zero-field (ZF) and transverse- 
field (TF) /iSR measurements in a single crystal of 
La2CuO4.11 (LCO:4.11) where excess oxygen is interca- 
lated in a stage-4 structure [16]. This system is supercon- 
ducting below T c <~ 42 K. Observation of sharp satellite 
Bragg peaks in neutron scattering indicates static long 
range (> 600 A) spin density wave (SDW) order below 
T N ~ 42 K. The modulation wave vector of the SDW is 
comparable to those observed in La2- a; Sr :c Cu04 (LSCO) 
and La 2 _ :E _ !y Nd 2/ Sr :I .Cu04 (LNSCO) systems with the 
hole concentration x ~ 0.125 [9,12]. The LCO:4.11 sys- 
tem is especially interesting since (a) it has the highest 
superconducting T c in the 214 family of materials; (b) it 
has a rather high magnetic TV, which is very close to T c , 
and well developed long range magnetic order below Tn] 
(c) there is no randomness arising from Sr substitution, 
and the intercalated oxygen ions are three dimensionally 
ordered. Preliminary results and analysis of our ZF-/uSR 
in LCO:4.11 have been presented at a recent conference 
[26]. 

We also report yuSR results in a superconducting 
Lai.88Sr .i2CuO4 (LSCO:0.12) single crystal with T c ~ 
30 K which exhibits magnetic Bragg peaks in neutron 
scattering [13]. Comparing these results with those of a 
simulation of the local field distribution at possible muon 
sites, we examine various models for the magnetically or- 
dered regions. The interplay between superconductivity 
and magnetism is discussed in the context of the super- 
fluid density n s /m* derived from TF-/zSR results. 

fiSR studies of magnetism in HTSC systems started 
in 1987 with studies of antiferromagnctic La2Cu04 (AF- 
LCO) [27], which provided the first evidence of static 
magnetic order. Subsequently, /xSR studies on La2Cu04 
[28], (La,Li) 2 Cu0 4 [29], La 2 _ x Sr x Cu0 4 [30-34], and 
YBa2Cu30j, [34,35], were performed to characterize mag- 
netic phase diagrams as a function of hole concentration 
and also to study details of the spin glass states near 
x = 0.05. 

/iSR measurements of HTSC systems with hole 
concentration x near 1/8 were first performed in 
Lai.875Bao.i 2 5Cu04 [36]. These detected static mag- 
netic order below = 32 K and a Bessel function 
/iSR line shape, characteristic of spin density wave sys- 
tems. A reduction of T c , associated with an increase 
of the muon spin relaxation due to quasi-static mag- 
netism, was also found by /iSR in La2- a; Sr a ;Cu04 near 
x = 0.12 [37]. Recent /iSR studies on (La,Nd,Sr) 2 Cu0 4 



(LNSCO) [38,39], (La,Eu,Sr) 2 Cu0 4 (LESCO) [40,41] 
and Lai. 8 7 5 Bao.i 2 5_ y Sr y Cu04 (LBSCO) [42] found the 
characteristic Bessel function line shape with the same 
frequency v <~ 3.5 MHz in all these systems. In the lat- 
ter studies the primary emphasis was determining the 
magnetic phase diagram as a function of rare-earth and 
hole concentrations. Although the existence of "zero/low 
field muon sites" , expected for decomposition of the sys- 
tem into regions with and without static magnetic order, 
was discussed in these reports [40,42], a systematic study 
of magnetic volume fraction was not made. 

Previously, Pomjakushin et al. [43] reported a /itSR 
study for single crystal specimens of La2Cu04+j, doped 
with oxygen, with y = 0.02 (LCO:4.02) and 0.04 
(LCO:4.04), which have compositions in the miscibility 
gap between AF-LCO and the stage-6 superconducting 
compound. In zero field, they observed /iSR precession 
signals identical to those found in AF-LCO. With de- 
creasing temperature in LCO:4.02, the amplitude of this 
oscillatory signal increased from nearly zero above the 
superconducting T c ~ 15 K to more than a half of that 
in AF-LCO at T ->• 0. In LCO:4.04, the oscillatory sig- 
nal appeared below T ~ 230 K with an amplitude corre- 
sponding to approximately a 10 % volume fraction for the 
AF region, and then exhibited a sharp increase below the 
superconducting T c ~ 25 K to about a half volume frac- 
tion at T — > 0. In both cases, the amplitudes of the AF 
oscillation exhibited an increase below T c . The diamag- 
netic susceptibility of the LCO:4.02 specimen was, how- 
ever, destroyed by a small applied field, indicating rather 
fragile superconductivity. Lack of information about the 
crystal orientation in the /ttSR measurements prevented 
a reliable estimate of the volume fraction in this study. 

In Section II we present experimental details and re- 
sults of the ZF-^SR measurements. Section III contains 
our simulation based on a model in which the sample 
contains microscopic regions where there is full magnetic 
order and other microscopic regions where the magnetic 
fields are too small to cause muon precession. In Section 
IV we presents results of measurements in a transverse 
magnetic field that are sensitive to the superconductiv- 
ity. Finally in Section V we discuss the results and draw 
conclusions. 



II. ZERO FIELD ^SR MEASUREMENTS 

Single crystals of LCO:4.11 and LSCO:0.12 have been 
prepared as described in refs. [16] and [13], which re- 
port neutron scattering studies of these crystals, respec- 
tively. The excess-oxygen doped La2Cu04+y sample is 
prepared by electrochcmically doping a single crystal of 
pure La2Cu04 which is grown by the traveling solvent 
floating-zone method. The crystal has a mass of 4.21 
grams and is cylindrical in shape. As result of twinning, 
there are equal populations of two twin domains with ei- 
ther the a or b crystallographic axis (in orthorhombic no- 
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tation) nearly parallel to the cylinder's long-axis. In both 
twin domains, the c-axis is perpendicular to the long- 
axis. Thcrmogravitimetric analysis on our samples has 
been performed on two smaller L^CnO^y crystals with 
T c = 42K. We find oxygen concentrations of y = 0.10(1) 
and y = 0.12(1). The large crystal has the same T c , 
so its oxygen concentration is expected to correspond to 
y = 0.11(1). _ 

Neutron diffraction measurements have shown that 
this material has a structural modulation along the c-axis 
corresponding to stage-4. [16]. Further details regard- 
ing the characterization of the stage-4 crystal by neutron 
scattering, susceptibility, and transport techniques may 
be found in references [16] and [44]. 

^iSR measurements were performed at TRIUMF, us- 
ing a surface muon beam with an incident muon mo- 
mentum of 28 MeV/c, following the standard procedure 
described in refs. [24,25]. We employed the so-called 
"low-background" muon spectrometer, which has the ca- 
pability of vetoing events from muons landing in areas 
other than the specimen, to ensure much less than ~ 10 
% background signal in our p,SR data. In ZF-^iSR, we ob- 
served positron time spectra via counters placed forward 
(F) and backward (B) to the incident beam direction (z), 
with the polarization of incident muons P M anti-parallel 
to z. The muon spin relaxation function G z (t) was ob- 
tained from forward (F) and backward (B) time spectra 
as 

G z (t) = [B(t)-F(t)}/[F(t)+B(t)\ (1), 

after correcting for the difference in effective solid angles 
of the F and B counters. As illustrated in Fig. 1(a), 
two sets of measurements were made with the c-axis of 
the specimen mounted parallel (configuration CZ-I) and 
perpendicular (CZ-II) to z \\ P^. 

Figures 2(a) and (b) show the muon spin polar- 
ization function G z (t) observed in these two crys- 
tals for CZ-I and CZ-II. Also included are the 
results obtained for ceramic specimens of non- 
superconducting Lai.875Bao.i25Cu04 (LBCO:0.125) [36] 
and Lai.475Ndo.4Sro.i3Cu0 4 (LNSCO:0.13) [38]. 

We see almost no relaxation above TV in the time 
range t < 1 (is. Below T N ~ 42 K of LCO:4.11 and 
Tjy ~ 20 K of LSCO:0.12, a small oscillatory component 
appears, with a Bessel-function line shape characteristic 
of muon precession in SDW systems [45,46]. The am- 
plitude A or of this oscillating and relaxing signal is lim- 
ited to about 34 % (CZ-I) and 24 % (CZ-II) at T 
for LCO:4.11, and 20 % (CZ-I) and 10 % (CZ-II) for 
LSCO:0.12. If we define the internal field at a muon site 
as H and 9 is the angle between z and H, as illustrated in 
Fig. 1(b), the oscillatory / relaxing amplitude represents 
sin 2 (#) averaged over all the muon sites. In a ceramic 
specimen which undergoes magnetic order in 100 % of 
the volume, this amplitude is 67 % (=2/3), as shown 
in the case of LBCO:0.125. Therefore, the observed sig- 
nal amplitudes A or for LCO:4.11 and LSCO:0.12 are less 



than a half the value expected for a sample with a SDW 
uniformly established throughout its volume. 

In single-crystal specimens, A or can have a small value 
if the local fields H for the majority of muon sites are 
parallel to z || P^. For this case, however, one should 
see a large A or when the crystal orientation is rotated by 
90 degrees from CZ-I to CZ-II or vice versa. Thus, our 
results, with small A or values for the both configurations, 
cannot be attributed to such anisotropy. In view of the 
symmetry of the LCO and LSCO structures, the "powder 
average" value of A or - pw d can be obtained [26] as 

-^or—pwd = 

(1/3) x [A or: cz-i + 2A or ..cz-ii], (2), 

which yields 27 % for LCO:4.11 and 12 % for LSCO:0.12. 
The reduction of A or ^ pwc i from 66.7 % indicates the ex- 
istence of "zero/low-field" muon sites. 

We have analysed the ZF-/xSR data by fitting the re- 
sults at t < 1/is to 

G z {t) = A or j (wt)exp(-At) + A n exp[-A 2 t 2 /2] (3), 

where j a is the 0-th order Bessel function and A n repre- 
sents the non-oscillating amplitude. The first exponen- 
tial term represents the damping of the Bessel oscillation, 
whereas the second describes the slow decay due to nu- 
clear dipolar fields. We obtain the volume fraction of 
muon sites that experience the SDW field using 

V„ = A or _ pwd /(2/3), (4) 

with A or _ pw d determined from Eq. (2) at the lowest tem- 
perature. We then scale the values for higher tempera- 
tures using A or results for the CZ-I configuration. For 
our fit range of t < 1/is, those muon sites with H > 30 
G contribute to the first term in Eq. 3, and thus arc 
identified as "finite-field muon sites". 

Figures 3(a) and (b) show the temperature de- 
pendence of and the oscillation frequency v = 
u/2tt for LCO:4.11, LSCO:0.12, LBCO:0.125 [36] and 
LNSCO:0.13 (ceramic) [38]. The common Bessel- 
function line shape and the nearly identical frequency 
v ~ 3.5 MHz (corresponding to maximum internal field 
of 260 G) at T — > for all the four systems indicate 
that the spin configurations and magnitudes of the static 
Cu moments (proportional to v) surrounding the "finite- 
field muon sites" are identical for all the four systems. 
The volume fraction increases progressively below 
Tjv with decreasing temperature, in both LCO:4.11 and 
LSCO:0.12. The frequency v acquires its full value just 
below Tjy in LCO:4.11. In Fig. 3(c), we compare the 
temperature dependence of the neutron Bragg-pcak in- 
tensity I B of LCO:4.11 [16] with V» x v 2 from /zSR. The 
good agreement of neutron and muon results indicate 
that these two probes detect the same static magnetism 
even though the former measures the ordered moment 
over distances > 600 A. The /iSR results suggest that 
the observed temperature dependence of Ib should be 
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ascribed mainly to the change of the site fraction con- 
taining ordered Cu spins, rather than to the increase of 
the static moment on individual Cu atoms although this 
is model dependent. 

To illustrate further this unusual temperature depen- 
dence of the magnetic order parameter in LCO:4.11, 
we compare ZF-yuSR spectra of LCO:4.11 and AF-LCO 
(T N > 250 K) in Fig. 4. In AF-LCO (Fig. 4(b)), as 
in most conventional magnetic systems where long-range 
magnetic order develops in 100 % of the volume, ZF-/zSR 
data exhibit muon spin precession whose amplitude is 
nearly independent of temperature, but whose frequency 
increases with decreasing temperature. In contrast, for 
LCO:4.11 (Fig. 4(a)) the oscillation amplitude increases 
gradually with decreasing temperature below Tn while 
the frequency just below Tn is already close to the low- 
temperature value. This implies that in LCO:4.11 the 
local order is well developed but dynamic above Tn, and 
the ordering process below T N involves the ordering of 
this local magnetism over long distances with a concomi- 
tant freezing of the spins. 

III. SIMULATION OF ZERO-FIELD ^SR 

To explore the possible origin of the two kinds of muon 
sites, we have performed computer simulations of the lo- 
cal field distributions at possible muon sites. The loca- 
tion of the muon site in AF-LCO has been approximated 
by Hitti et al. [47] to be about 1 A away from the apical 
oxygen, based on the expected length of the 2 ~-^ + hy- 
drogen bond. We have improved their estimate by taking 
into account the tilt of the CuOg octahedra and by com- 
paring the observed /iSR precession frequency v = 5.8 
MHz with the Cu long-range ordered moment size 0.55 
/i b in AF-LCO [16]. Figure 5 shows the location of the 
apical-oxygen muon site at (0.199<i, Od, 0.171c) in the 
unit cell with d = 3.779 A, and c = 13.2 A assumed in 
this simulation. Since the direction of H is very sensitive 
to a small change of the tilt angle of the CuOg octahe- 
dra, we calculate the expected ZF-fiSK line shape for the 
powder-average. 

Figure 6(a) shows the line shape for this apical-oxygen 
muon site for: (1) a 100 % volume fraction of antifcrro- 
magnetically (AF) ordered Cu moments having modula- 
tion vector kAF with ordered moment of 0.55 [Ib (chain- 
dotted line); and (2) Cu spins with an incommensurate 
SDW amplitude modulation superimposed on the anti- 
ferromagnetic correlations, with a maximum static Cu 
moment value of 0.35 hb (broken line). Here we have as- 
sumed that the SDW propagation vector is along the Cu- 
O-Cu bond direction (i.e., 45 degrees to the orthorhombic 
a and b axes) with magnitude k = 2ir/d x 0.12 (d = 3.78 
A is the distance between nearest neighbour Cu atoms) 
We have calculated the spin polarization using 

cos(kAFx) x cos(fcaf) = 



(1/2) cos[{k AF + k)x] + (1/2) cos[{k AF - k)x] , (5) 

at a given position x. These results confirm that the 
SDW model gives a line shape nearly identical to a Besscl 
function (solid line). In AF-LCO the observed /xSR re- 
sults [27,47] show a substantial damping of the sinusoidal 
precession, due to various possible origins, such as nu- 
clear dipolar fields, variations of the muon site or crystal 
imperfection. To account for the reduction of the ob- 
served ^iSR frequency from 5.8 MHz in AF-LCO to 3.64 
MHz at T — > in LCO:4.11, we need to assume that the 
maximum static Cu moment size in the SDW modulation 
is ~ 0.36 /i b ■ This may be compared with the bare Cu 
moment of 1.1 /xg. 

Since there is a substantial density of intercalated oxy- 
gen in LCO:4.11, we have also modeled the field distribu- 
tion at several locations 1 A away from the expected site 
of the intercalated oxygen atoms. As may be seen from 
Fig. 6(b), which displays decay curves calculated for 100 
% ordered SDW Cu moments with maximum moment 
size of 0.35 /ib, there is no muon site location having 
zero or low internal field. Therefore, we cannot ascribe 
the existence of "zero/low- field sites" in LCO:4.11 to in- 
tercalated oxygen. Furthermore, in LSCO:0.12, there is, 
of course, no intercalated oxygen, so it seems clear that 
the non-magnetic sites must have a different origin. 

In Fig. 7(a), we show the line shape G z {t) calculated 
for Cu moments in circular islands of SDW order in the 
Cu02 planes, for various volume fractions Vc u of static 
Cu spins, assuming the radius R of the circular islands 
to be R ~ 50 A. Bessel function line shapes are ob- 
tained with a common frequency, while the regions with- 
out static Cu moments account for the zero-field sites. 
The size of the magnetic island is reflected in the damp- 
ening of the Bessel oscillation. 

In Fig. 7(b) we compare the results of simulations 
for two different geometries. One is for circular magnetic 
islands with V Cu = 50% and R = 15 and 100 A. The other 
corresponds to an "ordered sandwich" model, where two 
out of every four Cu02 planes contain completely ordered 
SDW Cu spins while the other two planes, presumably 
those adjacent to the intercalated O layers, have no static 
Cu spins. The curve for R — 15 A islands shows a fast 
dampening of the Bessel oscillation, together with a slow 
decay of the non-oscillating signal. In contrast, the curves 
for R = 100 A islands and for the sandwich model exhibit 
a long-lived Bessel oscillation, with little decay of the 
non-oscillating signal. 

After generating line shapes for the SDW-island model, 
we have fitted the simulation results to the functional 
form of Eq. 3. Figure 8(a) shows the resulting mag- 
netic site fraction as a function of volume fraction 
containing static Cu spin order Vc u - We see that in- 
creases linearly with increasing Cu fraction, with a slope 
that is independent of island size R. The site fraction 
Vfj, becomes nearly 100 % at Vc u ~ 75%. This implies 
that even the observation of a magnetic /xSR signal cor- 
responding to — 100 % cannot rule out the existence 
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of regions of non-ordered Cu moments in ~ 25 % of the 
volume. By using the relationship between V M and Vcu, 
we estimate as lower limits Vcu ~ 27 % for our crystal 
of LCO:4.11 and - 10 % for LSCO:0.12. Of course,' this 
same possible difference between and Vc u applies in 
all other copper oxide SDW systems. 

In Fig. 8(b), we show the dampening rate A of the 
Bessel function oscillation. We find that a smaller radius 
R of the island results in a larger A. We also fit the 
observed spectra in LCO:4.11 with Eq. 3 and plot the 
corresponding point in Fig. 8(b) (open star symbol) . The 
result A4.11 = 1.25/is _1 agrees well with the simulation 
for R = 15 A. In actual systems, however, additional 
factors cause dampening of the Bessel oscillation, such 
as variations in the muon site, nuclear dipolar fields, and 
the effects of imperfections. 

To account for these additional effects, we have fit the 
spectrum of AF-LCO at T = 20 K in Fig. 4(b) with the 
simple form cos(wi + (f>) multiplied by an exponential de- 
cay exp(— At), and obtain Aaf ~ 0.81/zs . Assuming 
that Aap represents the additional relaxation contribu- 
tions, we have obtained a corrected relaxation rate from 

Al 11 = (Ai 11 -Ai F )V2. ( 6 ) 

That is, we assume that the decay constants add in 
quadrature. 

A4 n is shown by the closed star symbol in Fig. 8(b). 
These considerations lead us to estimate R = 15 ~ 30 A 
for the average radius of the static magnetic islands in 
LCO:4.11. A similar analysis for LSCO:0.12 would give 
roughly the same estimate for R. However, because the 
oscillation amplitude in LSCO:0.12 is much smaller than 
that in LCO:4.11, the analysis would be less reliable. 

It is also possible to obtain the relaxation rate expected 
for the sandwich model. The result (A = 0.2/xs -1 ), 
shown in Fig. 8(b) by the asterisk symbol, indicates 
that this model is less successful than the magnetic is- 
land model in reproducing the relaxation A4 n observed 
in LCO:4.11. 

We now discuss a simple heuristic model for the propa- 
gation of magnetic order between these presumed islands. 
If we assume that spins in neighbouring islands in the 
same CuC>2 plane are correlated when there is any over- 
lap in the islands' illustrated in Fig. 9, we can 
estimate the effective correlation length of the spin order. 
For various values of Vc u we have calculated the proba- 
bility that the direction of the Cu spins in the magnetic 
islands, located at a distance r from a Cu spin in the cen- 
ter of a magnetic island, is correlated with the direction 
of the Cu spin at r = 0. The results of this calculation, 
shown in Fig. 9, suggest that even with the assumption 
that the spins are in islands, the order in a given plane 
may nevertheless propagate over long distances. How- 
ever, for Vcu = 30% the correlation length from Fig. 9 
is only ~ 3 times the size of the magnetic island, whereas 
the narrow SDW peaks observed in neutron scattering 
measurements indicate that the static spins are corre- 
lated over very large distances (> 600 A) within the plane 



[16]. Therefore, some additional mechanism to increase 
the connectivity is required. 

The results in Fig. 9 are obtained assuming no spin 
correlations between magnetic islands belonging to dif- 
ferent neighbouring Cu02 planes. The neutron measure- 
ments of LCO:4.11 [16] indicate that the static SDW or- 
der exhibits short range correlations along the c-axis di- 
rection with a correlation length of ~ 13 A, somewhat 
larger than twice the distance between adjacent Cu02 
planes. Such spin correlations along the c-axis direction 
would further increase the size of the correlated clusters 
shown in Fig. 9. For example, if we allow correlations 
when any overlap exists among areas of the islands pro- 
jected to the adjacent island on neighbouring planes, the 
effective number of islands contributing to the "connec- 
tivity" would increase at least by 3 times compared to 
the case without interplane correlations. This is because 
all the islands on a given Cu02 plane plus those on the 
upper and lower planes would participate in the correla- 
tions. 

In Fig. 10, we illustrate this by showing randomly po- 
sitioned islands of 30 % integrated area fraction in three 
different planes (a)-(c). When we overlap these planes, 
almost all the islands achieve percolation, as shown in 
Fig. 10(d). Moreover, there could be further contribu- 
tion from islands on the second neighbour planes which 
would enhance connectivity within the first neighbour 
planes. In this way, we expect a strong tendency to- 
wards percolation of spin correlations among randomly 
positioned magnetic islands. 

IV. SUPERCONDUCTING PROPERTIES 

TF-^SR measurements allow us to derive the mag- 
netic field penetration depth A from the superconduct- 
ing contribution a of the muon spin relaxation rate as 
a oc A~ 2 oc n s /m* (superconducting carrier density / 
effective mass). This rate a reflects the inhomogeneity 
of the magnetic field in the flux vortex structure, which 
varies over a length scale of ~ 1000-3000 A. Any het- 
erogeneity in the superconducting properties at length 
scales shorter than this will be averaged out. Thus, /zSR 
probes superconductivity with a ~ 100 times coarser spa- 
tial resolution than that with which it probes static mag- 
netism. In HTSC systems without any static magnetic 
order, the line shapes in TF-^iSR spectra can be used 
to study more detailed spatial features, such as the size 
of the vortex core region [48]. In the present case, with 
co-existing superconductivity and static magnetic order, 
however, detailed line-shape analyses become very diffi- 
cult. 

We have performed TF-/xSR measurements by rotating 
the incident muon polarization to be perpendicular to the 
beam P M _L z and applying an external field parallel to z 
with the c-axis of the specimen mounted either parallel 
to z [configuration CT-I] or perpendicular to z [configu- 
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ration CT-II], as illustrated in Fig. 11. The muon spin 
polarization function G x {t) is measured using two sets 
of counters placed up (U) and down (D) of the incident 
beam. 

As shown in Figs. 12(a) and (b), we observe a larger 
relaxation in CT-I, the geometry reflecting the in-plane 
penetration depth A a &, than in CT-II. We have analyzed 
the TF-^iSR results by assuming the existence of a static 
internal field in a volume fraction of muon sites, hav- 
ing the distribution in magnitude and direction estimated 
from the ZF-^SR measurements. We have then added 
vectorially the external field, with a Gaussian broaden- 
ing due to superconductivity. For the "zero/low- field 
muon sites" with volume fraction (1-V^), we assume that 
the muon sees only the nuclear dipolar fields and the 
Gaussian-broadened external field. We fit the observed 
data with this model, and derive the Gaussian half-width 
at half maximum of the superconducting contribution 
equal to cr/j^. Figure 12(a) shows the calculated con- 
tributions to the relaxation from superconductivity and 
static magnetism, separately. 

The temperature dependence of a, observed in the CT- 
I configuration, is shown in Fig. 13 for LCO:4.11 and 
LSCO:0.12. The relaxation rate shows a gradual increase 
with decreasing temperature below T c . To compare the 
values of a(T — > 0) with those observed in ceramic spec- 
imens of other HTSC systems [49-52], we multiply by 
1/1.4 to account for the effect of the anisotropic penetra- 
tion depth [53] . We then add the corresponding points to 
the plot of a{T — > 0) versus T c in Fig. 14. The points for 
the present systems lie on the same trajectory as other 
LSCO systems. We have also included a point for LESCO 
(EuO.l, Sr0.15; ceramic), which exhibits static magnetic 
order with ~ 50% [40,54,55]. This point again falls 
on the trajectory. The correlation in Fig. 14 suggests 
that a oc n s /m* at T — > is a determining factor for 
T c in HTSC systems with static stripe freezing (present 
systems and LESCO) in a way similar to the case for 
underdoped [49,50], Zn-doped [52] and overdoped [51,56] 
HTSC systems. 

V. DISCUSSION AND CONCLUSIONS 

A. Interplay between magnetic and superconducting 
volumes 

Comparing the results of a to those observed in ce- 
ramic LSCO systems, assuming that m* is independent 
of doping, the results for LCO:4.11 and LSCO:0.12 cor- 
respond to volume average hole densities of ~ 0.14 +/- 
0.02 and 0.10 +/- 0.02 holes per Cu, respectively. In view 
of the large systematic errors, however, we can not use 
these results to distinguish whether the volume without 
static magnetism (1-Vcu) alone or the total volume car- 
ries superconductivity. We can rule out, however, a case 
where superconducting carriers exist only in the mag- 



netic volume Vcui since the local hole-concentration for 
this model would need to be unrealistically large (more 
than 0.3 holes per Cu) to account for the observed value 
of volume integrated n s /m* 7 which is comparable to 
those in LSCO systems in the optimum-doping region. 
So far, a clear indication of mutually exclusive regions 
with static magnetism and superconductivity has been 
obtained only in ^iSR measurements in (La,Eu,Sr) 2 Cu04 
(LESCO) [40,54,55], which demonstrate that the super- 
fluid density n s /m* scales as (1-V^). 

In previous work on Lai.4 5 Ndo.4Sr .i5Cu04 and 
Lai.45Ndo.4Sro.2Cu04 [38], which are superconducting 
below T c ~ 7 K and 12 K, respectively, ~ 100 % was 
found. This result could be interpreted as evidence of 
spatial overlap of regions supporting superconductivity 
and regions having static spin order. The relationship 
in Fig. 8(a), however, indicates that Vc u f° r this sys- 
tem can be smaller than 100 %. This leaves open the 
possibility that superconductivity with reduced T c sur- 
vives in a small remaining volume fraction (l-Vcu) ~ 20 
%, while superconducting regions and regions with static 
magnetism are mutually exclusive. 

The situation with magnetic islands, having length 
scales comparable to the in-plane superconducting coher- 
ence length, resembles the "swiss cheese model" [52] for 
Zn-doped HTSC systems, where a non-superconducting 
region of comparable size is created around each Zn. This 
model was proposed based on the ^SR results for the re- 
duction of n s /m* as a function of Zn concentration [52], 
and was recently confirmed by direct measurements of 
scanning tunnelling microscopy (STM) [57]. The /iSR re- 
sults for overdoped HTSC systems can also be explained 
if one assumes spontaneous formation of hole-rich non- 
superconducting islands embedded into a sea of hole-poor 
supcrfluid [55, 58-60]. The correlation between T c and 
n s /m* survives in all these systems, as shown in Fig. 14. 
We note that similar microscopic heterogeneity in the su- 
perconducting state has also been found in recent STM 
measurements on underdoped Bi2212 systems [61]. 

The T c vs. n s /m* correlations are robust against var- 
ious perturbations in HTSC materials, such as those 
caused by Zn-impurities, overdoped fermion carriers, and 
the formation of static magnetic islands, seen in the 
present work. This is analogous to the robustness of cor- 
relations between the superfluid transition temperature 
and the 2-dimensional supcrfluid density in thin films of 
4 He and 3 He/ 4 He adsorbed in porous and non-porous 
media [62-66]. Based on these observations, one of us 
[55,59,60] pointed out a possible relevance of these het- 
erogeneous electronic/magnetic features in HTSC sys- 
tems to a "microscopic phase separation", such as the 
one seen in superfluid 3 He/ 4 He films adsorbed in porous 
media / fine powders [66,67]. 

In the present work on LCO:4.11, as well as a previous 
,uSR study on LCO:4.02 and LCO4.04 [43], the onset of 
superconductivity occurs around the temperature below 
which the volume fraction of static magnetism increases 
with decreasing temperature. This apparent coincidence 
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of T c and TV could be expected if a microscopic rear- 
rangement of charge distribution occurs at T c ~ TV to 
separate the electron systems into hole-rich regions which 
support superconductivity and hole-poor (x ~ 1/8 for 
LCO:4.11 and x ~ for LCO:4.02 and LCO:4.04) re- 
gions which support static magnetism. 

Recently, Kivelson et al. [68] have proposed a model 
with a particular type of phase separation to explain the 
sharp rise of v and gradual increase of Vc u below Tjv 
seen in the present study of LCO:4.11. Their model also 
predicts a trade-off of superconducting and magnetic vol- 
umes below T c ~ Tjv- The gradual change of both <r(T) 
and Vv(T) below T c - T N in LCO:4.11 (Figs. 3(a) and 
13), however, seems inconsistent with such a trade-off. 

The superconducting T c and the magnetic TV ap- 
pear at different temperatures in uSR measurements of 
LSCO:0.12 as well as in various LNSCO and LESCO sys- 
tems, with LSCO:0.12 having T c > TV and the other ma- 
terials having T N > T c . It should be noted that neutron 
measurements give T c ~ Tn in LSCO:0.12, while T/v de- 
termined by uSR is lower, presumably due to difference 
in time windows between neutron and muon measure- 
ments. In these systems the frequency v increases with 
decreasing temperature gradually below T/v, as shown in 
Fig. 3(b), and in refs. [38-42]. Thus, the abrupt develop- 
ment of the magnetic order parameter is not a common 
feature of all the HTSC systems. 

B. Bragg peak intensity in LCO:4.11 and AF-LCO 

We can estimate the Bragg peak intensity Iba.u ex- 
pected in neutron scattering measurements on LCO:4.11 
from the /xSR results. Iba.u should be reduced from 
the value Ib-.af in AF-LCO by a factor 0.5 for the SDW 
amplitude (sin 2 ) times the ratio of the maximum static 
moments (0.36/0.55) 2 times the volume fraction of mag- 
netic Cu atoms Vc u = 0.27. Thus, we predict as a rough 
estimate Ib-a.ii = 0.06Ib-.af- This value agrees well with 
the observed neutron results Iba.u = 0-07Ib-.af- In Ref- 
erence [16] it is assumed that three out of four copper 
ions are magnetic and that they are distributed uniformly 
throughout the entire volume, so a smaller average mo- 
ment per Cu 2+ of 0.15 \ib is inferred. 

We note, however, that estimates of the ordered mo- 
ment size from neutron studies and from uSR often dis- 
agree with one another. For example, the Neel tem- 
perature in antiferromagnetic I^CuCU+a varies from 
T/v ~ 300 K to Tn < 100 K with a small variation in 
S. In neutron scattering studies of these AF-LCO sys- 
tems, the Bragg peak intensity Tb shows a reduction, by 
more than a factor of 10, with decreasing T/v [28]. In 
contrast, the frequency of ZF-/xSR spectra shows only a 
20 % change in the local frozen moment for the same 
reduction of T/v- These /xSR and neutron results on AF- 
LCO are compared in ref. [28] . Since the /xSR frequency 
is directly proportional to the magnitude of neighbouring 



static Cu moments, the reduction of Is in neutron studies 
cannot be ascribed to a change of individual Cu moment 
size. The likely reason for the strong reduction of Is is a 
progressive trade-off between the true long range ordered 
component and short range spin glass fluctuations which 
recent work [69] has shown involves fluctuations into the 
diagonal stripe spin glass phase with hole concentration 
of -0.02. 



C. Neutron results in high magnetic fields 

Recently Khaykovich et al. [48] have made measure- 
ments on LCO:4.11 in high magnetic fields H ext applied 
parallel to the c-axis. By increasing the external held 
from H ext = to H ext = 9 T, they have found a factor 
— 2 increase of the intensity Tb of the magnetic Bragg 
peak at T — ► for a crystal similar to the one stud- 
ied here. The superconducting T c is suppressed by H ext , 
while Tn remains nearly unchanged up to H ext = 9 T. 

Application of H ext perpendicular to the Cu02 plane 
creates vortices. The vortex core, with radius comparable 
to the in-plane coherence length £ a b — 30 A becomes 
normal, as illustrated in Fig. 15. These normal cores 
could have static magnetic order similar to that in the 
surrounding magnetic islands. 

The fraction of the area contained in vortex cores may 
be estimated from the upper critical field, since that is 
the field at which the cores fill the entire area. Using 
H C 2 > 40 T, found for optimally doped LSCO [70], we 
expect that only — 1/4 or less of the area that is super- 
conducting at H ext — would be turned into a normal 
core at H ext = 9 T. However, this is sufficient to explain 
the factor — 2 increase in the Bragg peak intensity if the 
zero-field sample has ordered moments in only ~ 27% of 
its volume and if the magnetism in the cores is coherent 
with the long range order present at zero field. We are 
currently undertaking /iSR studies of LCO:4.11 in high 
applied magnetic fields, which will be reported in future 
publications. 

In nearly optimally doped LSCO, Lake et al. [71] have 
recently found that low-energy neutron scattering inten- 
sity, below the energy transfer corresponding to the su- 
perconducting energy gap, increases at temperatures well 
below T c , when a high external field H ext is applied along 
the c-axis. This phenomenon may also be related to the 
magnetic response from the flux-core regions, since the 
increase of the intensity sets in at the flux-depinning tem- 
perature. A rather surprising feature of this study is the 
sharpness of the observed low-energy fluctuations in mo- 
mentum space, suggesting long-range spin correlations 
among widely separated vortex core regions. We point 
out here that, like the static order, this requires long 
range coupling between vortex cores. 
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D. Conclusions 

In summary, we have performed ZF-^iSR measure- 
ments using LCO:4.11 and LSCO:0.12 single crystals, 
and have found that static incommensurate SDW spin 
freezing, with the maximum ordered Cu moment size of 
0.36 /is, develops only in a partial site fraction. We as- 
sume Vf, ~ 40% for LCO:4.11 and - 18% for LSCO:0.12, 
and specific modeling suggests that the corresponding or- 
dered Cu fractions may be somewhat smaller. Compar- 
ison of observed results with computer simulation sug- 
gests the formation of static magnetic islands on the 
CuC>2 planes having size R — 15 ~ 30 A, comparable to 
the in-plane superconducting coherence length £ a b. Or- 
der between these islands nevertheless propagates over 
long distances to yield a correlation length in excess of 
600 A. 

We stress that LCO:4.11, which is a stoichiometric sin- 
gle crystal, has the least built-in randomness among the 
various HTSC systems that exhibit static incommensu- 
rate magnetic correlations. Yet the present work demon- 
strates that the ground state is a mixture of magneti- 
cally ordered and at best weakly ordered regions. Near 
phase boundaries or quantum critical points in strongly 
correlated electron systems, such microscopic heterogene- 
ity could result from competition of the different order 
parameters involved. Related phenomena have been ob- 
served, for example, in the formation of stripe magnetic 
correlations in manganites [72] and the development of 
stripe domains in LSCO with x ~ 0.01 to 0.02 [73]. 
These observations encourage further theoretical studies 
of electronic states involving spontaneous formation of 
heterogenous regions in competing order parameter sys- 
tems. 

Our TF-^iSR measurements in these systems demon- 
strate that n s /m* at T — * exhibits correlations with T c 
not only in underdoped, Zn-doped, and overdoped HTSC 
systems but also in the present systems with static SDW 
spin freezing. Recent measurements of the penetration 
depth in (BEDT-TTF) 2 Cu(NCS) 2 in applied pressure 
[74], and in A 3 Ceo systems [75,76], suggest that these 
correlations are followed also by organic and fullerine su- 
perconductors. 

We have proposed a model for explaining the rather 
long-range spin correlations in these systems, resorting 
to "connectivity" of neighbouring magnetic islands on 
the same and adjacent Cu02 planes. We have also sug- 
gested that static magnetism in the vortex core region 
can provide explanations for the field dependence of the 
neutron scattering results in LCO:4.11 and in optimally 
doped LSCO systems. 

The present work demonstrates a very important fea- 
ture of /iSR, i.e., the capability to determine the volume 
fraction of magnetically ordered regions. Recently, a sim- 
ilar case was noticed in the study of URu 2 Si2, where 
static magnetism below TV ~ 17 K has been identi- 
fied to exist in a partial volume fraction by combination 



of neutron [77] and NMR [78] studies in applied pres- 
sure. The first evidence for this feature was indeed pro- 
vided by a jLtSR measurement [79] in ambient pressure. 
In the ^SR results for URu 2 Si2, Luke et al. found a 
very sharp onset of the precession frequency v(T) be- 
low Tjv, temperature dependence of the precessing am- 
plitude U M (T), and a good agreement between the tem- 
perature dependence Ib(T) of the neutron Bragg peak 
and VpV from /iSR. These features are common to the 
present results in LCO:4.11. We also note that ZF-/uSR 
results in CeCu2.2Si2 [80,81] indicate a possible trade- 
off between magnetic and superconducting volume frac- 
tions. These results suggest possible involvement of mi- 
croscopic phase separation in both HTSC systems and 
heavy fermion systems [55,68,82]. Further detailed stud- 
ies of ,uSR in combination with neutron scattering will 
be very helpful in obtaining an overall understanding of 
the interplay between magnetism and superconductivity 
in HTSC, heavy-fermion, and other strongly correlated 
electron systems. 
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FIG. 2. Time spectra of the muon spin polarization 
observed by ZF-/iSR with (a) the CZ-I configuration, and 
(b) the CZ-II configuration for La2CuO4.11 (LCO:4.11) and 
Lai.88Sro.i2Cu04 (LSCO:0.12). Also included are results for 
ceramic specimens of Lai. 875 Bao.i25Cu04 (LBCO:0.125) from 
[36] and Lai.47Ndo.4Sro.i 3 Cu0 4 (LNSCO:0.13) from [38] 



FIG. 1. (a) Schematic view of the experimental configu- 
rations in zero field ^iSR (ZF-^SR) employed in the present 
study. Measurements were performed with the initial muon 
spin polarization parallel (CZ-I) and perpendicular (CZ-II) to 
the c - axis, by rotating the crystal orientation by 90 degrees. 
The upper [lower] figure shows the CZ-I [CZ-II] configuration, 
(b) Illustration of the oscillating part of the muon polariza- 
tion signal, proportional to sin 2 (8), with 8 being the angle 
between the local magnetic field and the initial muon spin di- 
rection. The cos 2 (8) component contributes to the amplitude 
of the non-oscillating signal. 
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FIG. 3. (a) Volume fraction V M of muon sites with 
a static magnetic field larger than ~ 30 G and (b) fre- 
quency of the precessing signal in La2CuO4.11 (LCO:4.11) and 
Lai.88Sro.i2Cu04 (LSCO:0.12) (from the present study), com- 
pared with the results in Lai.s75Bao.i25Cu04 (LBCO:0.125) 
[36], Lai.47Ndo.4Sro.i3Cu0 4 (LNSCO:0.13) [38], and antifer- 
romagnetic La 2 Cu0 4+ 5 (AF-LCO) [27,28]. The broken lines 
are guides to the eye. (c) Comparison of the temperature de- 
pendence of the Bragg peak intensity Ib of neutron scattering 
measurements in La2CuO4.11 (LCO:4.11) [16] with those ex- 
pected from the /u,SR results (present study) as Ib oc V m x v 2 . 
/iSR and neutron results are scaled using the values near 
T — > 0. Comparison of the absolute values is discussed in 
Section V-B. 
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FIG. 4. Time spectra of ZF-/iSR measurements in (a) 
LCO:4.11 and (b) AF-LCO, at various temperatures. The 
amplitude of the oscillating and relaxing signal varies in (a) 
below Tjv without much change in the frequency, while in (b) 
the amplitude does not depend on T and the frequency in- 
creases with decreasing T. The latter behavior is observed in 
ZF-/1SR of many conventional magnetic systems. The spec- 
trum (a) was obtained by using low-background /iSR spec- 
trometer, which ensures that the background signal level is 
below 0.01 in Asymmetry. The spectrum (b) was obtained in 
1987 by using older apparatus. The broken line (BG) indi- 
cates an expected level of background signal in (b). 
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FIG. 5. Schematic view of the crystal structure of 
La2Cu04, shown with the positions of intercalated interstitial 
oxygen and the muon site associated with an apical oxygen. 
We assume the location of muon site near apical oxygen to be 
(0,0.751,2.257) A within the unit cell. The arrows attached 
to O 2- in the center of the figure show the direction of tilting 
(about 5°) of the CuC>6 octahedra. 
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FIG. 6. Computer simulation results of the muon spin 
polarization functions expected for La2Cu04 assuming that 
all the Cu moments involved in the static magnetic order. An 
angular average is taken to represent the results for ceramic 
specimens, (a) Results for the muon site near the apical oxy- 
gen (see Fig. 5) . The dotted line corresponds to the case with 
antiferromagnetic (AFM) spin correlations with the static Cu 
moment of 0.55 /is, while the dashed line corresponds to the 
case with spin density wave (SDW) amplitude modulation 
with maximum static Cu moment of 0.35 hb and modula- 
tion wavevector k = 0.12 * 27r/d (d=3.779 A) superimposed 
on the antiferromagnetic correlations. The solid line is the 
result of a fit to a Bessel function, (b) Results for several 
hypothetical muon sites near the intercalated interstitial oxy- 
gen, illustrated in the inset, obtained for the above-mentioned 
SDW spin correlations. 
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FIG. 7. Computer simulation results of the expected 
muon spin polarization functions for La2Cu04, with only a 
part of the Cu moments involved in the static magnetic order, 
obtained for the muon site near the apical oxygen after an- 
gular average to simulate the results for ceramic specimens. 
We assume a SDW modulation with a maximum static Cu 
moment of 0.35 /is and modulation vector k = 0.12 x 2n/d. 
(a) Magnetic island model where the Cu moments in a vol- 
ume fraction Vcu order with the SDW amplitude modulation, 
forming islands of radius R — 50 A. (b) Comparison of the 
results for the magnetic island model with R —15 and 100 A 
with those for the "sandwich model" where two CuC>2 planes 
adjacent to the intercalated oxygen layer do not have any 
static Cu moments while all the Cu moments on the other 
two planes are ordered. All three curves shown in (b) corre- 
spond to Vcu=50%. The insets illustrate these two models, 
with the shaded regions containing SDW modulated static Cu 
moments and the blank regions without static Cu moments. 
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FIG. 8. (a) Computer simulation results for the volume 
fraction V M of muon sites with a static magnetic field larger 
than ~ 30 G as a function of the volume fraction of the sample 
containing static Cu moments, Vcu- (b) The relaxation rate 
A of the Bessel function oscillation, obtained by fitting the 
simulation results with Eq. 3, plotted as a function of Vcu- 
Comparison with the experimental results for LCO:4.11 (Raw 
data A4.11 shown by an open star symbol and corrected data 
A-4.il by a closed star symbol) allows estimation of the size of 
magnetic islands. The asterisk symbol * shows the relaxation 
rate A expected for the sandwich model. 
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FIG. 9. Probability P(r) that Cu spins belong to a certain 
cluster versus the distance from the center of an of radius R in 
the cluster. We assume that all the spins in a given island are 
correlated, and the correlation extends if there is any small 
overlap between adjacent islands on a given CUO2 plane, as 
illustrated in the upper figure. For example, a Cu spin at a 
distance r 2 contributes "1.0" for the probability, while that at 
r\ gives "0.0" contribution. We average over all the possible 
cluster/island configurations to calculate P(r). The results 
demonstrate a long length scale of spin correlations for the 
"magnetic island model" , due to formation of large clustered 
islands. 
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FIG. 10. Illustration of percolating cluster islands, (a)-(c) 
show planes with random locations of magnetic islands having 
integrated area fraction of 30 %. (d) shows the overlap of 
(a)-(c). (d) demonstrates that most of the islands belong 
to the percolating cluster if we allow correlations of spins in 
"overlapping islands" on all the three planes. 
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FIG. 11. Schematic view of experimental configurations 
of Transverse Field /iSR (TF-^iSR) measurements employed 
in the present work with the external field applied parallel 
(CT-I) and perpendicular (CT-II) to the c-axis of the single 
crystal. 



f— SC and SDW 
--SC only 




LC0:4.11 H =1kG 1 c T=3K 

ext 

I I I I 



400 600 

Time (ns) 



1000 



FIG. 12. Time spectra of the muon spin polarization ob- 
served in TF-/xSR in LCO:4.11 obtained using (a) CT-I and 
(b) CT-II configurations. The relaxation of the signal is due 
to a distribution of magnetic fields inside a superconducting 
sample (SC) and magnetic fields due to static Cu moments 
with the SDW amplitude modulation, (a) shows faster damp- 
ing than (b), reflecting the shorter penetration depth for the 
external field applied perpendicular to the Cu02 planes. The 
expected effect from superconductivity alone is shown by the 
dotted line in (a). 
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FIG. 13. Temperature dependence of the superconduct- 
ing relaxation rate a due to superconducting flux vortices, 
observed in TF-/1SR measurements (with the CT-I config- 
uration) in La2CuO4.11 (LCO:4.11) and Lai.8sSro.i2Cu04 
(LSCO:0.12). The broken lines are guides to the eye. 
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FIG. 14. A plot of the superconducting transition tem- 
perature T c versus the relaxation rate o(T — » 0) at low tem- 
peratures (proportional to the superfluid density n s /m*) for 
several high-T c superconductors [49-52]. The results with 
the "stripe square" symbols represent points from LESCO 
([54,55]), LSCO:0.12, and LCO:4.11 in the order of increas- 
ing a(T — > 0). To account for difference between results 
for ceramic and single-crystal specimens, the values for a for 
LCO:4.11 and LSCO:0.12 in the CT-I configuration were mul- 
tiplied by 1/1.4. 
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FIG. 15. An illustration of magnetic islands (shaded cir- 
cles) and vortex cores (open circles) on the CUO2 planes 
of HTSC systems under an external magnetic field applied 
along the c axis direction. When a core is located at the 
magnetic island, one expects an increase in the radius of 
the non-superconducting region due to supercurrent flowing 
around such a magnetic island. 
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